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For a class of nonnegative, range-1 pair potentials in one dimensional continuous space we prove that 
any classical ground state of lower density > 1 is a tower- lattice, i.e., a lattice formed by towers of 
particles the heights of which can differ only by one, and the lattice constant is 1. The potential may 
be flat or may have a cusp at the origin; it can be continuous, but its derivative has a jump at 1. The 
result is valid on finite intervals or rings of integer length and on the whole line. 
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q ! 1 Introduction 



A curious aspect of soft potentials is that they allow classical particles to superimpose in ground state 
configurations (GSCs), even if the pair potential is purely repulsive pQ. This property is related to 
the Fourier transform of the interaction and can appear — mostly but not exclusively — if this one is 
partially negative. Indeed, recently it was proven that a strictly positive Fourier transform leads to a 
uniform distribution of particles at high densities, while a partially negative Fourier transform causes 
non-uniformity 2 . However, the forms of non-uniformity can be various, and the precise conditions on 
the Fourier transform giving rise to superimposed particles are still unknown. The three-dimensional 
example of Likos et al. p] is the pair potential e~ a \ x \ with m > 2. In the case m — 4 a detailed 
numerical analysis and analytic arguments show rather convincingly the superposition of particles. The 
second derivatives of these interactions vanish at the origin. This implies a strong negative part of the 
Fourier transform at large wave vectors, and intuition also suggests that the flatness of the interaction 
at zero distance favors the formation of towers of particles at lattice sites. The type of the lattice is 
less intuitive. If the (negative) minimum of the Fourier transform decides the lattice type, the reciprocal 
lattice should be close-packed. Until now a rigorous proof exists only in the case of the penetrable-sphere 
model, particles with a repulsive square core potential [2]. This is the flattest possible interaction, its 
Fourier transform is oscillating and slowly decaying, and not the reciprocal but the real-space lattice 
is close-packed, the lattice constant being the range of the potential. Mathematically, this problem is 
almost trivial and, because of the jump of the interaction, the model is pathological in the sense that the 
energy and the free energy do not tend to the energy of the GSCs as the temperature goes to zero [2]. 
Below we present a one-dimensional generalization of the penetrable-sphere model, providing less trivial 
examples. We shall also see that a vanishing second derivative of the interaction at zero distance is not 
a prerequisite for the phenomenon to occur (this was noted already in [1]). The extreme value minus 
infinity for the second derivative, i.e., a cusp at the origin is also allowed. This is particularly interesting 
in comparison with other examples of interactions which have a cusp at zero and a partially negative 
Fourier transform, and no Bravais lattice, with or without superposition of particles, as a high-density 
GSC 0. 



2 Theorem and proof 

Given an even translation invariant pair potential u, the potential energy of N pointlike particles in the 
configuration (x)n '■= (xi,X2, ■ ■ ■ , %n) is 

U(x) N = ^ u(xi-xj). (1) 

l<i<j<N 

Further on, {x} m will denote the sequence x, . . . , x where x is repeated m times. Recall that an iV-particle 
GSC in a bounded domain A is a configuration that minimizes U(x)n among (x)n G A w , and a GSC 
in inhnite space is an infinite configuration whose energy cannot be decreased by any number-preserving 
local perturbation [2]. 

Theorem 1 In one dimension, let u be a bounded real function with 

( = 1, x = 

u(x)=u(-x)( >l-|x|, < \x\ < 1 (2) 

1=0, \x\ > 1. 

(i) GSCs on an interval of integer length. 

1. For arbitrary positive integers n, m, any r 6 {0, 1, . . . , n} and any integers < i\ < ■ ■ ■ < i r < n 
the configurations 



l—lr 



({hr + \ . . . , {i r r + \ !{./!"■*, ,) (3) 



are N = m(n + 1) + r-particle GSCs of u on the interval [0, n]. That is, each integer point in [0, n] 
is occupied by m orm + 1 particles. The energy of the N -particle GSCs is 



E n {N) = (n + 1) Q + rm = U(X^J. 



(4) 



2. Ifr = 0, the unique N -particle GSC in [0,n] is X n - m = ({j} m )™ =0 . 

3. If u{x) > 1 — \x\ for < \x\ < 1 then A™'™ are the only N-particle GSCs in [0,n] for every r. 

4- If u is the 'overlap potential', u{x) = (1 — |a;|)lr_i .i] (x), and r > then the ground state is contin- 
uously degenerate. On the background of X n,m the remaining r particles can be freely distributed 
under the constraint that the distance of any two of them is > 1 . 

(ii) GSCs on the whole line. 

1. Any infinite configuration X — ({i} mi )^° = _ 00 with rrii > 1 and \rrii — rrij\ < 1 for all i, j is a GSC 
of u on R. 

2. The ground state of the overlap potential is continuously degenerate. For any integer m > 1, any 
< Xi < . . . < x m < 1 and any Y = (• • • < y_i < yo < y\ < ■ ■ •) such that yk+i — Vk > 1 for all k, 



m oo 



X=\J |J {i + Xj} and XUY (5) 

j — 1 i— — oo 

are GSCs. Here the union is in the sense of superposition, i.e. coinciding points are multiply 
counted. 

3. If a GSC has a density, and the value of this is p, then the ground-state energy per unit length is 

lim ^ 1 = \[P\{P + { P }-1) (6) 

ri-i-oo n Z 

where \p\ and |p) are the inteqer and fractional parts of p, respectively. 



(Hi) GSCs on a ring of integer length. 

For n > 1 let A — [0,n+ 1[ and u\{x) — J2j = _ ao u(x + j(n + 1)). Any N > n + l-particle GSC of u 
on [0,n] is a GSC of u\ on A with energy E n (N). In the case of the overlap potential the restriction to 
[0, n + 1[ of any GSC of the type f5|) is also a GSC of u\. 

Proof, (i/1). Let < x\ < X2 < ■ • • < xn < n be any iV-particle configuration: 

U(x) N = E u(xj - Xj) > u(xj - x i+ j). (7) 

l<i<j<N i>l j=l 

If (x) jv = X™' i then (JT]) holds with equality because Xi + j — xi > 1 for j > m + 1. The number of 
retained pairs is 

/m + 1\ 

> 1,1 < j < m,i +j < N} = N - 1 + N - 2 H \- N - m = mN - . (8) 
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The lower bound (JT]) can be rearranged into a summation over chains of pairs, 



m m J 

EE"^ - x l+j ) = EE N 3 ^ - x i+ j) + u(x i+j - Xi+2j) + u(x i+2j - x i+3j ) H ]. (9) 

i>l J — 1 j — 1 z— 1 

There are altogether 53j=iJ = ('"J 1 ) chains. Now we apply the lower bound 

u(x — y) > 1 + x — y, 

valid for any pair x < y, to each term of ((9]). The sum of the unities gives the total number of pairs 
([5]). The sum over the differences of particle coordinates in a chain gives the smallest minus the largest 
coordinate, and can be bounded from below by —n. Thus, 

(m + 1\ (m + 1\ . , „ , . , (m + 1\ 
U(x)n > mN — ( 2 )-«( 2 J = m[m(n + 1) + r] - (n+ 1) f g 

(10) 



(n + l)^)+rm = f/(X^ ir ,) 



where the last equality is obvious. This shows that X^ m i is a GSC and proves the formula (U]) for the 
ground-state energy. 

(i/2). Let r = 0, i.e. TV = m(n + 1). If X n ' m = (£) N , & < then 

m— 1 

Lf(r">EE M te-^) ( n ) 

and, instead of ([5]), it suffices to sum over j up to m - 1. In 

>l,l<j <m-l,i + j <N} 
the number of pairs is (m — 1)^ — (™) and the number of chains is (™). This yields the lower bound 

Y^u(x t -x %+] )>{m-l)N-{n + l)(™^ = (m-l)m(n+l)-(n+l)M = (n+l)(^j = E n (N). 

(12) 

We show that if (x)n ^ X n,m , then there is a pair k + m) such that x/c+ m — x& < 1, and therefore 

m — 1 

c/(*v > E E 

) + w(:r fc - x k+m ) > E n (N) + 1 + x k - x k+m > E n (N). (13) 



Assume that Xk+ m — x% > 1 for all k. Then 

n 

n > sXxjm+k ~ X(i-l)m+fc) = x nm+k - X k > n (14) 
i=l 

for k — 1, 2, . . . , to. Thus, 

Xi = #2 = • • ■ = x m = and x nm +i — x nm +2 = • • ■ = X( n +i)m = n - 
It follows that x m+ i > 1 and x„ m < n — 1. Then 

n-l 

77 — 2 > / ^(Xj TO -|-fc — ^(i-l)m+fe) — a; (n-l)m+fc — ^m+fc > 77 — 2 (15) 
i=2 

for fc = 1, 2, . . . , to. Thus, 

Xm+l — X m +2 = ■ ■ ■ = X2m = 1 an d 3J(n— l)m+l = X(n-l)m+2 = ■ ■ ■ = X nm = Tl — 1 

and it follows that X2m+i — 2 and x< n -i)m < 77 — 2. Further repeating this argument, in \n/2\ steps we 
find that the unique (x)n m [0, n] satisfying Xk+ m — Xk > 1 for all k is X n,m . 

(i/3). Consider the case when u(x) > 1 — \x\ for < \x\ < 1. We may suppose r > 0. Let (x)jv be 
such that < Xfe — a;^ < 1 for some pair i < k. If k > i + to, the positive contribution of ti(xj — Xk) was 
not taken into account in (J7J). If fc < i + m, m(x^ — x^) was taken into account with the lower bound 
1 + Xi — Xk- In either case we obtain U(x)n > E n (N). Thus, in any GSC, for any i < k we have either 
Xk — X{ > 1 or Xi — Xk with k < i + to. Necessarily, Xk+ m +i — Xk > 1 for all k. As in (i/2), one can then 
prove that Aj"'™ are the only configurations satisfying this property. 

(i/4). If u(x) = 1 — | a; | for |a;| < 1 and is zero otherwise then 

U({i} m ,{x},{i + l} m ) = m 2 , (16) 

whenever x G + 1], and the particle in a; does not interact with particles not in i or i + 1. This proves 
the assertion. 

(ii/1). Take any A = ({i} mi )°° = _ ocl with to^ > 1 and Ito^ — raj| < 1 for all i,j. Thus, can assume 
only two values, say, m(> 1) and m + 1. Any local modification of X is confined in an interval [ri\, 712] 
where rti < n 2 are integers. According to (i), An [ni,n 2 ] is a GSC in [rii,n 2 ], meaning that for any 
N = \X n [ni, 71,2] (-particle configuration (y)./v C [711,712], 

£%V>£/(An [774,77,2]). (17) 

Because the distance of both (y) n and An [ni , 772] to the rest of A is at least 1, there is no interaction 
between particles inside [771,772] and in A \ [771,772]- Together with (|17[) this implies that A is locally 
stable, that is, a GSC onK|. 

(ii/2). This could be derived from Eq. (|16p. but there is a more transparent geometric proof: u{x) = 
(1 — \x\ )1 [-1,1] (x) is called the overlap potential because it measures the overlap of two aligned rods of 
unit length such that their centers are at a distance |x|. Replace each particle of A or of A U Y by a rod 
(interval) of unit length centered at that particle. Then no local displacement of the rods can decrease 
the sum of their overlaps, showing that A and A U Y are locally stable. 

(ii/3). For sufficiently fast decaying interactions the energy density, if exists, is the same (and minimal) 
for all locally stable configurations of the same density. This absence of metastability was proven in [3] . 
We can therefore compute the energy density from any GSC of a given density. Choose A = ({i} mi )^° = _ oa 
with \rrii — rrij\ < 1. Let / and £(I) denote a finite interval and its length, respectively. The lower and 
upper densities of A are defined as 

p(X) = liminf J^i^l, p(X)=limsup^L^ (18) 
where the limit is taken on arbitrary sequences of intervals. Suppose that 



o(X) = p(X) = pe [to, 777 + If. 



Let / = [0, n] and N = \X H [0, n] | . If n is large enough then N = m(n + 1) + r where < r < n and r 
depends on /. However, because of the existence of the limit defining the density, the limit of E n {N)/n 
also exists and equals ([5]). 

(iii) This is the case of a periodic boundary condition, and ua is the periodized interaction. Because 
the range of u is 1, in ua only a single term can be nonzero. The difference and relative subtlety compared 
with case (i) comes from the fact that now the length of the interval is larger by 1 for the same number 
of particles, but the interaction between particles at the two ends of the interval adds to the energy. Here 
we refer to a more general result. Fix N > n + 1 and let Xq be one of the TV-particle GSCs of u on [0, n]. 
Let X be the periodic extension of Xq to K. By part (ii), X is a GSC of u on R. Then, according to 
Lemma 7.2 of Ref. ||, Id A=Xq is a GSC of ua on A. Furthermore, 



(x,y)CX 

In the case of the overlap potential we can directly restrict the configurations ([5]) to A and find the result 
for example by the geometric argument of (ii/2). 

3 Discussion 

The aim of this paper was to provide provable nontrivial examples of pair interactions which give rise 
to the superposition of particles in classical ground state configurations. Obviously, such an interaction 
must be bounded for touching particles. Trivial examples include unstable interactions, for examples, 
attractive ones. An almost trivial example is the penetrable sphere model. The important discovery of 
Likos et al. [T] was that the interaction can be repulsive everywhere outside the origin, what really counts 
is that it must have a partially negative Fourier transform. This led them to the study of the family 
e -«N w ith m > 2. Although the superposition has been established numerically, a mathematical proof 
is still missing for this class of potentials. On the other hand, in [5] counterexamples were constructed, 
showing that the partial negativity of the Fourier transform was not sufficient — it was, however, proven 
to be 'almost necessary' in the sense that a strictly positive Fourier transform causes the particles to 
distribute uniformly as the density increases. The examples given in the present paper are nonnegative 
range- 1 potentials which may or may not be purely repulsive within their range. The GSCs at densities 
p < 1 are trivial and form a continuously degenerate family (any configuration in which the distance 
of neighboring particles is larger than or equal to 1 is a GSC), while in GSCs for p > 1 the particles 
superimpose on the sites of a lattice of lattice constant 1. It is quite possible, however, that for 'generic' 
interactions having a partially negative Fourier transform the high-density GSCs are not tower-lattices 
but periodic configurations in which the particles cluster around (and not on) the sites of a lattice. 

Potentials whose Fourier transform is nonnegative and takes on zero somewhere represent a marginal 
case. They can be obtained as limits of potentials having a partially negative Fourier transform. If these 
latter have tower-lattice GSCs, the property extends to the limit by continuity, but the tower-lattice is 
expected to be embedded in a continuum of other GSCs. Here, this is the case of the overlap potential 
(previously studied by Torquato and Stillinger [4 ) whose Fourier transform is (2//c 2 )(l — cos k). At integer 
densities the ground state of the overlap potential on an interval of integer length is non-degenerate, but 
in infinite space and on the ring we meet the same kind of continuous degeneracy as in the case of pair 
potentials with a nonnegative Fourier transform of compact support [5]. 

There has been no use of Fourier transform in our proof. This may seem curious in view of the 
decisive role of the Fourier transform of the interaction. However, in fc-space a proof would have been 
less complete. Based on Proposition 4.2 of Ref. [2J, for integer densities p = m (and only for them) we 
can prove that X n ' m are GSCs of % on A = [0,n + 1[ by showing that among positive measures p of 
total weight n + 1 on A the functional 





(19) 




n+1 



n+1 




(20) 



is minimized by po — Sj=o $o (here Sj is the Dirac delta localized at j). In Fourier representation I[p] 
reads 




where 



1 



I 



n+l 



u(k) — / u(x)e 



—ikx 



(22) 



n + l 



Therefore, one must prove that 



oo 



u(0) _ 1 



- 6(2i7r) 



2 



- 2 <m- 



(23) 



For the overlap potential I\po] is easily seen to be the minimum, and because the other potentials of the 
family @ are bounded below by the overlap potential, I[^o] is the minimum of I[u] also for them. 
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